Leggett mode in a strong-coupling model of iron arsenide superconductors 



F. J. Burnell, 1 ' 2 ' 3 Jiangping Hu, 4 Meera M. Parish, 5 ' 6 and B. Andrei Bernevig 1 ' 5 

Department of Physics, Princeton University, Princeton, NJ 08544 
2 Kavli Institute for Theoretical Physics, Santa Barbara, CA 
3 All Souls College, Oxford, OX14AL 
^Department of Physics, Purdue University, West Lafayette, Indiana 47907 
5 Princeton Center for Theoretical Science, Princeton University, Princeton, NJ 08544 
6 Cavendish Laboratory, J J Thomson Avenue, Cambridge, CBS OHE, UK 
(Dated: April 22, 2010) 

Using a two-orbital model of the superconducting phase of the pnictides, we compute the spectrum 
of the Leggett mode - a collective excitation of the phase of the superconducting gap known to 
exist in multi-gap superconductors - for different possible symmetries of the superconducting order 
parameter. Specifically, we identify the small regions of parameter space where the Leggett mode 
lies below the two-particle continuum, and hence should be visible as a sharp resonance peak. We 
discuss the possible utility of the Leggett mode in distinguishing different momentum dependencies 
of the superconducting gap. We argue that the observation of a sharp Leggett mode would be 
consistent with the presence of strong electron-electron correlations in iron-based superconductors. 
We also emphasize the importance of the orbital character of the Leggett mode, which can result in 
an experimental observation of the mode in channels other than A± g . 



I. INTRODUCTION 

The discovery of high-temperature superconductivity 
in iron arsenide and related compounds at the begin- 
ning of 200& 1 has triggered an enormous interest in the 
condensed matter physics community and has stimulated 
a flurry of experimental activity^! Upon electron^ 
or hole 7 doping of a magnetically-ordered parent state, 
most of the iron-based superconductors exhibit transition 
temperatures T c beyond the conventional BCS regime, 
with some extending up to 56 K 11 , thereby breaking the 
cuprate monopoly on high-temperature superconductiv- 
ity. Experimental evidence accompanied by theoretical 
modeling suggest that the pairing in the iron-pnictides is 
different from the d-wave pairing of the cuprates. Never- 
theless, they resemble the cuprates in that it is increas- 
ingly clear that the magnetism of the parent state (either 
long-range or fluctuating order) crucially influences the 
pairing symmetry of the doped system. A conclusive ob- 
servation of the pairing symmetry still remains elusive, 
with both nodal and nodeless order parameters reported 
in experiments. This provides a strong incentive to iden- 
tify new experimental probes potentially sensitive to the 
symmetry of the superconducting gap. 

While a wide r ang e of nodal gap functions were ini- 
tially predicte d 12 " 14 !, the general theoretical view has 
now converged to favor an extended s-wave order pa- 
rameter (denoted ,s ± or s x 2 y 2) that takes opposite signs 
on the electron and hole pockets along the multi-band 
Fermi surfaces. The symmetry of this s x 2 y i gap matches 
that of the iron-pnictide Fermi surface: it is maximal 
around (0, 0), (7r, 0), (0, 7r), (7r, tt) - the location of the 
Fermi surfaces in the unfolded one-iron-per-site Bril- 
loiun zone. This sign- alternating nodeless gap is con- 
sistent with some experimental data and also has b road 
theo retical sup portF^l] Indeed, both stron^pDIL anc i 
weaM^^- 6 * 18 * 20 " ^- coupling theories of the onset of super- 



conductivity predict an extended s-wave order parame- 
ter. 

Experimentally, however, there is no consensus about 
the nature of the order parameter, with both nodal and 
nodeless gaps being reported. While most experiments 
can be explained within the framework of an s ± gappU, 
several facts, such as the T 3 dependence of the NMR re- 
laxation rate over a significant temperature range^MD, 
residual finite quasiparticle terms in the thermal conduc- 
tivity^ISl, as well as the power-law behavior of the pene- 
tration depth 33 ! 34 ! , remain unsettled. Some of the exper- 
iments on penetration depth and thermal conductivity 
could be explained by an s ± order parameter if there 
were a large gap anisotrop ^ 331341 , but this contradicts 
ARPES data, which reveals very isotropic nodeless gaps 
on the hole Fermi surfaces^^, of magnitudes matching 
a strong-coupling form A(fc) = Ao cos(k x ) • cos(k y }^ in 
the unfolded Brillouin zone. 

A possible resolution of this apparent contradiction, 
consistent with the theoretical prediction of an s± or- 
der parameter, is that the gap anisotropy is doping de- 
pendent and that different experiments are done at dif- 
ferent dopings. In the strong-coupling mean-field pic- 
tur d 17 * 19 * , the gap anisotropy is intrinsically doping de- 
pendent: the gap has a form cos(k x ) • cos(k y ) which be- 
comes more anisotropic as the doping is increased. In a 
weak-coupling expansion of Fermi surface interactions, 
the gap anisotropy can arise from the presence of an 
Ai g term cos(k x ) + cos(k y ) (which does not break the 
crystal symmetry but can create nodes on the (it, 0) and 
(0, 7r) electron surfaces) in the band interactions 3 -^ upon 
renormalizatiori 2 ^. A large gap anisotropy is already 
present in Fu nctio nal Renormalization Group studies of 
orbital models^HEll 

In this paper, we analyze another physical phenomenon 
- the Leggett mode of multi-band superconductors - that 
depends on the strength of the pairing order parameter 
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and could also in principle quantitatively distinguish a 
sign- changing gap from other gap symmetries. Specifi- 
cally, we investigate to what extent the pairing symmetry 
of the iron-based superconductors can be deduced by an- 
alyzing the behavior of the Leggett mode as a function of 
doping and the strength of superconducting order param- 
eters. As the iron-based superconductors have multiple 
orbitals, the superconducting state exhibits a plethora of 
collective modes beyond the usual Goldstone/Higgs plas- 
mon. Here we use an effective two-orbital model of the 
superconducting state to study one of these - the Leggett 
mode associated with anti-symmetric phase fluctuations 
between the two superconducting order parameters. This 
gapped collective mode can, in the right parameter range, 
present a sharp collective mode resonance below the two- 
particle continuum which could in principle be detected 
experimentally. 

To determine whether such a collective mode reso- 
nance occurs in the pnictide superconductors, we study 
the gap and dispersion of the Leggett mode as a func- 
tion of doping and the superconducting order parame- 
ters. We show that, for a sign-changing gap function, the 
Leggett mode can be below the two-particle continuum 
for a small regime at low doping. In particular, when 
the band renormalization is large, an undamped Leggett 
mode can exist in a relatively large parameter region. 
Thus, the observation of a sharp Leggett mode will vali- 
date the presence of strong electron-electron correlations 
in the iron-based superconductors. Moreover, in our two- 
orbital model, the Leggett mode is a B\ Q mode, instead 
of a pure A\ g mode, which is expected in any band-based 
model. Therefore, the orbital structure of pairing in the 
iron-based superconductors can be validated by identify- 
ing the existence of the Leggett mode in channels other 
than Ai g . 

Unfortunately, we find that the Leggett mode cannot 
qualitatively distinguish between a sign-changing order 
parameter and other gapped order parameters. However, 
the sign-changing order parameter will have a degree of 
anisotropy which depends on doping. For large doping, 
the sign-changing order parameter on the enlarged Fermi 
surface will exhibit larger anisotropy. As such, the super- 
conducting gap will be small on some parts of the Fermi 
surface and the Leggett mode will be overdamped, lying 
above the two-particle continuum, and hence unobserv- 
able. This presents a testable opportunity if, at moderate 
doping (when the gaps should theoretically be isotropic) , 
the Leggett mode is below the two-particle continuum 
and hence observable. If so, the observation of a disap- 
pearing collective mode provides indirect support for a 
sign-changing gap function. 



II. COMPUTING THE EFFECTIVE ACTION 

The Leggett mode is a collective excitation of two- 
(or multi-) band superconductors, associated with anti- 
symmetric phase oscillations between the two bands. It is 



thus a neutral mode associated with oscillations between 
the supercurrents of the two bands. Here we present the 
effective action for this mode derived from a two-orbital 
model appropriate to the pnictides at temperatures well 
below the onset of superconductivity. To render our cal- 
culations analytically tractable, we focus on a simplified 
model of the iron-based superconductors that takes into 
account only the d XZl d yz orbitals. In this case, an intra- 
band order parameter can have its phase fluctuate be- 
tween the two orbitals in two modes: the usual sym- 
metric combination (Goldstone) and the antisymmetric 
combination, which is the Leggett mode. 

While the conventional Leggett mode involves only the 
Fermi surface gaps, our work involves a Leggett mode in 
the orbital gaps. We consider the orbital basis rather 
than an effective band basis, because neglecting the or- 
bital structure of the iron-based superconductors is most 
likely incorrect: it was showrP^that due to the difference 
in mirror symmetry eigenvalues of the electron and one of 
the hole bands at the Y point in the Brillouin zone (BZ), 
the spin density wave (SDW) state is gapless with a Dirac 
point in both two and five-orbital models of iron-based 
superconductors. This highly nontrivial effect, confirmed 
by experiment^, is lost in the effective band-basis pic- 
ture. Details of the derivation of the Leggett mode ef- 
fective action in the orbital basis, which differs slightly 
from the band-basis result of, e.g., Ref. [4TJ are given in 
Appendix [Aj 



A. Model Hamiltonian 

Using the insight provided by numerical and ana- 
lytic studies suggesting that the antiferromagnetic ex- 
change coupl ing between next-nearest-neighbor Fe sites 
is strong 12 42 , two of us^ studied a t-J\-J2 model with- 
out band renormalization and obtained a gap function of 
the form <zos(k x ) • cos(/^), which changes sign between 
the electron and hole-pockets of the Fermi surface of the 
material. It is this type of strong-coupling superconduc- 
tivity that we will focus on in this paper, but we point 
out that other weak coupling approach es exist an d give 
a similar sign-changing order paramete J 15 * 16 * 18 ^^. 

To calculate the effective action for the phase modes 
of the superconducting state, we employ a model of the 
pnictides which incorporates only the d xz and d yz or- 
bitals at each site, together with hybridization between 
the two. Although this description is only truly valid in 
the case of unphysically large crystal field splitting, we 
use this model for its analytic simplicity. We adopt the 
band structure proposed in Ref. [14[ which at first glance 
captures the essence of the Density Functional Theory 
results: 



where %jr k a = (c* d fc(T ,c^ k a ) is the creation operator 
for spin a electrons in the two orbitals and the kinetic 
terms read: 

£x(k) = —2t\ cos k x — 2t2 cos k y — 4ts cos k x cos k y — fi 
e y (k) = —2^2 cos k x — 2t\ cos k y — 4t 3 cos k x cos k y — fi 
txy(k) = — 4^4 sin k x sin k y (2) 

The hoppings have roughly the same magnitude: t\ = 
-1.0, t 2 = 1.3, t 3 = -0.85, and U = -0.85 in eV. We find 
that the half-filled, two electrons per site configuration is 
achieved when \i = 1.54eV. 

The missing ingredient in this two-orbital model is the 
d xy orbital, which can be shown to be important to the 
detailed physics of the iron-based superconductor^. For 
example, the kinetic model ([I]) gets the location of the 
second hole pocket wrong - it situates it at the (n, tt) 
point in the unfolded BZ, whereas LDA calculations show 
two hole pockets at the T point. However, the two-orbital 
model gets several of the qualitative characteristics of the 
iron-based superconductors right: it has a nodal SDW 
instability and a sign-changing s-wave superconducting 
instability. 

To describe the superconducting phase, we use the 
approach of Ref. 17, adopting a strong-coupling pic- 
ture in which the interaction Hamiltonian contains 
anti-ferromagnetic nearest-neighbor and next-nearest- 
neighbor coupling between the spins in both identical 
and opposite orbitals. While not entirely correct at lat- 
tice scales, it was shown that this model gives remarkably 
large overlaps with the interactions obtained through the 
functional renormalization groupP^ method, and hence 
can be considered as an effective interaction model for 
the iron-based superconductors. Furthermore, for our 
purposes these interactions are important only insofar as 
they give, after decoupling in the superconducting chan- 
nel, the sign-changing cos(k x ) • cos(k y ) superconducting 
order parameter. In this sense, the interacting spin- 
model we use can be thought of as an effective Ginzburg- 
Landau description of iron-based superconductors; the 
precise mechanism driving the transition to the super- 
conducting phase is irrelevant to the effective action we 
derive here. 

Based on the mean-field analysis of Ref. 17, we will as- 



sume throughout that the superconducting instability is 
dominated by the intra-orbital interactions, so that the 
gap is diagonal in the orbital basis. Indeed, at the mean- 
field level, the inter-orbital pairing is weaker than the 
intra-orbital pairing by a factor of approximately five^. 
In addition there is a large on-site inter-orbital Hund's 
rule coupling which will not enter into the present analy- 
sis as it does not alter the nature of the order parameter 
at mean-field level. We will briefly discuss the impacts of 
this last term, together with the antiferromagnetic inter- 
orbital interactions, in Sect. |II C| 

B. Phase-only effective action 

To obtain an effective action for the phase of the super- 
conducting gap, we follow the general protocol of Ref. l45| 
details of this calculation as applied to the orbital basis 
are given in Appendix [Aj In essence, one first decouples 
the interaction terms in the microscopic model using a 
Hubbard- Stratonovich transformation. This re-expresses 
operators quadratic in the fermions as interaction terms 
between a pair of fermions and the superconducting field 
<£. Deep in the superconducting region, where fluctua- 
tions in the magnitude A = |3>| can be neglected, inte- 
grating out the fermions then yields an effective action 
for the phase modes of the system. Since we work with a 
two-orbital model, there are a priori two superconducting 
gaps, excluding inter-orbital pairing. Though by symme- 
try their magnitudes have to be equal, this leads to two 
independent phase degrees of freedom. As is well known, 
one of these is a Goldstone mode which, upon including 
the Coulomb interactions, becomes a plasma mode. The 
other is the (gapped) Leggett mode, which will be our 
principle focus here. 

For our purposes, the two phase degrees of freedom are 
most conveniently expressed in the basis 

4,= -L(8 1 + 6 2 ) ^,= -^-02) ( 3 ) 

where 9\ and 62 are the phases of the gaps in the xz and 
yz orbitals, respectively. Hence <f> represents the symmet- 
ric phase oscillation, while <p represents the (neutral) an- 
tisymmetric phase mode. In this basis, we find the effec- 
tive action to be (see calculation details in Appendix [A}: 



S eJf = J dnd 2 q(<P(Sl,q) <p(Sl,q) 



n 2 

^2 



N vv [tf 



iQiQj 



n 2 



(4) 
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with momentum-independent coefficients given by: 



\2 



N U = - / 7^Zv2 { —7a^ + 737a^ f ( 5 ) 



? (A)3 

t 2 




(6) 



M = / — (7) 



with f2o = }J -rf ■ Here, E'i are the two band energies E± = ^ ( e x + e y ± ^/ (e x — e y ) 2 + 4e^ y j of the metallic state, 



and ££ A) = + A 2 are the quasi-particle energies in the superconducting phase. All e, E, and A are evaluated 

at the momentum k to be integrated over. The above equations represent the main result of the paper. 

In Eq. Q, terms linear in g, as well as terms bilinear in g, O, all vanish in the limit T — »• 0. As expected, this 
effective action Q describes one gapless mode, comprised entirely of symmetric phase fluctuations at q = 0, and one 
gapped mode. The latter is the Leggett mode; at q = it consists purely of antisymmetric phase oscillations between 
the two superconducting gaps. Here we are principally interested in the Leggett mode gap, ^o, as this represents the 
threshold at which the mode becomes experimentally observable. Thus if < 2A, we expect the Leggett mode to 
appear as a sharp resonance in the spectrum of the pnictide superconductors. 

For terms involving q 2 , the expressions for the coefficients c 2 a p - are somewhat more complicated and are thus given 
in Appendix Al We note, however, that for i ^ j, any coefficient of q^qj vanishes due to symmetry. Further, for 



i = j, symmetry of the coefficients under a 90 degree rotation of the Brillouin zone fully determines their direction 
dependence in q. Taking these symmetries into account, Eq. Q has the form: 



Q 2 _ r 2 a 2 




s eff = j xuP q ( ^ ) ^ J - - - y! > • (8) 



We should note that the Leggett mode gap is propor- 
tional to e 2 y - that is, to the off-diagonal kinetic terms 
in the orbital basis. This is in contrast to the approach 
of, for example, Ref. STJ in which the superconducting 
gap is taken to be diagonal in the band basis of the nor- 
mal state, and it is the inter-band interactions which 
couple the phases of the two gaps, and hence generate 
the Leggett mode. This difference stems from the fact 
that we take the gap to be diagonal in the orbital ba- 
sis: A a (k) = (c a ^kCai-k), where a indexes the orbitals 
and assume that the pairing is defined over the whole 
Brillouin Zone. Any model in which the interaction is 
written in orbital space and which aims to respect the 
point-group symmetries of the lattice will require this 
type of orbital-basis formalism. 

C. Including Hunds interactions 

In light of the fact that our approach is based on an 
absence of off-diagonal interactions in the superconduct- 
ing channel (when the orbital basis is used), it is useful to 
consider in more detail the validity of this assumption in 



the presence of inter-orbital couplings. In the pnictides 
the ferromagnetic Hund's rule interaction 

Hh = —Jh ^ Sl r S2r = —Jh ^ 4rcr a ^,^ /C lrcr / 4r-Y a 7,7 /C 2r 7 / 
r r 

(9) 

is the principal source of such interactions. 

Since spin ordering must be absent in the supercon- 
ducting phase, generically we may decouple the Hunds 
interaction in either the particle-particle channel or the 
particle-hole channel. At lowest loop order, the particle- 
particle interaction serves only to renormalize the band 
structure. The particle- hole contribution was, as previ- 
ously noted, shown to be small by Ref. [T7| Neglecting 
the small inter-orbital pairing at mean-field, we find that 
the Hunds interaction affects the effective action for the 
Leggett mode cp only through higher loop corrections in 
the fermion propagator. 

Further, it is straightforward to include the effect of 
the small inter-orbital interaction in the superconduct- 
ing channel. Such a term simply modifies the effective 
action for the superconducting phase by adding a term 
V12 (A1A2 + A 2 A^) = 2|A| 2 Vi2Cos((^). This modifies 



5 



the gap of the Leggett mode according to 

V 12 



V U V 2 



22 



V? 2 N vv 



(10) 



Here V a p parametrize the superconducting interaction 
between orbitals a and /?, as described in Eq. |A2[ 
and we have taken A(fc) = A cos{k x ) cos(k y ). For 
< V12 <C Vii, V22, the effect of including such a term is 
always to bring the Leggett mode gap down in energy. 



gle band superconductor, including the Coulomb inter- 
actions modifies the effective action for the phase 6 of 
the superconducting gap such that becomes a plasma 
mode 46 . We will not examine in detail the plasma mode 
here; rather we note that including Coulomb interactions 
does not substantially modify the relevant features of the 
Leggett mode, as we show below. This is not surprising 
since the Leggett mode is, at long wavelengths, associ- 
ated with the neutral antisymmetric phase oscillations, 
and hence does not couple to the Coulomb interaction. 



D. Effective action with Coulomb terms 

In the above analysis, we ignored the effects of the In the presence of Coulomb interactions, the phase- 
Coulomb interaction on the phase modes. In a sin- only effective action has the form: 



S eff 



J dVtd 2 k {(j) (p' 



l-U(q)N^ 




(11) 



(For the sake of completeness we derive this result in Appendix [bJ . The dispersion of the symmetric mode <\> is 
modified by the denominator of its Q 2 term; in practice since the Coulomb interaction U (q) is singular as q this 
makes the symmetric mode into a plasma mode, exactly as is known to occur in a single band super conduct oJ^. 

Eq. ( pTj ) shows that including Coulomb interactions does not alter the mass gap of the Leggett mode, as the plasma 
mode does not mix with the Leggett mode tp at q = 0. The net effect of the Coulomb terms on ip(£l, q) will be a 
modification of the q 2 term in the effective action of the Leggett mode. Integrating out </>, we obtain: 



eff 



dftd 2 k(p(q)(p(-q)N ip 



n 2 



n 2 



Ml 



Q 2 



A a V~lAp (12) 



Hence for small q, Q, in the presence of Coulomb interactions, provided that lim g ^ q^i-N^v^)) * s ^ m ^ e ' tne ne ^ 
effect is a modification of the effective velocity of the mode. The above equation needs to be solved self-consistent ly 
to obtain the mode dispersion. However, the limit q — >• 0, which determines whether the Leggett mode is above or 
below the particle-particle continuum, is unchanged from the case without the Coulomb interaction. 



III. RESULTS 



Having established the general form of the effective ac- 
tion of the Leggett mode, we now turn to a quantitative 
evaluation of the coefficients in Eq. Q. Our principle 
interest will be what potential information the Leggett 
mode can give about the form of the superconducting 
gap - in particular, we address the question of whether 
it can distinguish between the popular extended s-wave 
gap and other plausible pairing symmetries. Unfortu- 
nately, it is clear from our equations that the Leggett 
mode properties depend on the absolute value of the gap 
function, thereby preventing any qualitative sensitivity of 
the mode to a sign-change in the gap function. We find 
that the clearest signature is the lifetime of the Leggett 
mode as a function of doping - at low dopings we find the 



Leggett mode to lie below the two-particle continuum; at 
higher doping the mode is always at higher energies than 
the two-particle continuum and hence will give at best a 
very broad resonance. 

A. The Leggett Mode gap 

We begin by studying the Leggett mode gap f^o for 
several different gap functions, with the objective of un- 
derstanding the qualitative differences expected between 
these in potential experiments. In each case, the mode is 
expected to be visible if it lies below the two-particle con- 
tinuum, which is set by 2 min |A| (where the minimum is 
taken over the Brillouin zone). 

At q = and T = 0, the symmetric and antisymmetric 
phase oscillations decouple, and from the effective action 
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(a) (b) 




FIG. 1: Gap of the Leggett mode for (a) extended s-wave, (b) standard s-wave, (c)A = A sin fee sin/cy, and (d) cZ-wave 
superconducting gaps. The magnitude of Qo in eV is indicated by the color map to the right of each figure, with blue 
corresponding to regions of smaller f2o, and red to regions of larger Qo- The vertical axis indicates the filling zy, with 1/2-filling 
corresponding to the undoped case; the horizontal axis is the scale of the maximum magnitude of the gap in eV: we take 

A = A r fe . 



( 12 ) the gap of the Leggett mode cp is given by 



M 

1 * in in 



with M and given by Eqs. ^ and 0. Note that 
M > and iV w < 0, so that the Leggett mode gap is 
well denned. We can evaluate the coefficients and 
M by integrating the expressions ([5| and Q numerically 
over the Brillouin zone. We use the values of e a quoted 
in Eq. (p). 

Fig. 1 shows the expected gap of the Leggett mode 
for extended s-wave, standard s-wave, d-wave, and 
A sin k x sin k y gaps, as a function of the filling fraction 
v and the maximum gap magnitude Ao- The general 
form of is similar in all four cases: it increases with 
the superconducting gap Ao, and has its lowest values 
at a filling of approximately v = 0.4. For all four order 
parameters, we also find the gap of the Leggett mode 
shows an academically interesting chemical potential de- 
pendence, droping sharply between v = 0.4 and v = 0.5 



independent of the momentum-dependence of the order 
parameter. 

The qualitative features of these plots can be under- 
stood by considering the form of Eqs. ([5| and (|7|). First, 
we see that increases monotonically with Ao, at a 
slightly less than linear rate. Though naively both M 
and N w scale quadratically with Ao, N w has diver- 
gences if E^ or vanish; these are cut off by the 
gap but nevertheless contribute the major part of the 
integral. Consequently, N w is well approximated by 



Vfs'/A, with Vfs the volume of the Fermi sur- 



face. On the other hand, M vanishes at the Fermi surface 
in the limit of small A, scaling approximately as M A in 

this region. Hence the quantity 



M 



increases with 



A, with a power close to (but slightly less than) 1. 

The non-monotonic dependence on v, which is similar 
in all four cases, stems from the dependence of the shape 
and volume of the Fermi surface on the chemical poten- 
tial. As stated above, the integral expression for N w is 
dominated by contributions near the Fermi surface. M, 
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FIG. 2: Plots of Fermi surface as a function of chemical potential for 0.3 < v < 0.65. The relevant filling fractions are shown 
as v in the title of the figure. The electron pockets first appear at approximately v — 0.4, and cross the nodal lines at v — 0.6. 
The nature of the hole pockets does not change substantially over the range shown here. 



on the other hand, receives significant contributions from 
the entire Brillouin zone, and is thus much less sensitive 
to the shape of the Fermi surface. 

To illustrate how these effects play out in our two- 
orbital model, Fig. [2] plots the Fermi surface of the two- 
orbital model in the normal state for the range of dopings 
considered here. Between v = 0.4 and v = 0.45, at pre- 
cisely the locus of the sharp drop in Qq seen above, a new 
set of Fermi pockets appears at the points (±7r, 0) and 
(0, ±7r) as a second pair of bands crosses the Fermi level 
in the normal state. In the superconducting state this 
results in more areas where the expressions for are 
relatively large - in particular, due to the much smaller 
Fermi velocity near the new branches of the Fermi sur- 
face, the area over which N w is large increases sharply, 
leading to a sudden reduction in Qq. Also worthy of note 
are the extremal values (y = 0.3, ^ > 0.6) at which the 
Fermi surface intersects with the nodes of the extended s- 
wave gap. These account for the non- monotonic behav- 
ior of observed in both extended s-wave and d-wave 
order parameters between v — 0.3 and 0.4, as the cut- 
off in the normal-state divergences of grows smaller. 
Though the application of our simple two-orbital model 
at large fillings is not warranted, and the features dis- 
cussed in this paragraph are model-dependent, we expect 
them to be accurate for gaps diagonal in the orbital basis 
inasmuch as the band structure given by the two-orbital 
model is correct. 



B. Observability of the Leggett mode 

In order for the Leggett mode to give a sharp resonance 
in experiments, it should lie below the two-particle con- 
tinuum. For the <i-wave and sine- wave gaps, which are 
nodal for the iron-based superconductors' Fermi surfaces, 
this is obviously never the case. For ordinary s-wave and 



extended s-wave gaps, the position of the Leggett mode 
at q = relative to the two-particle continuum depends 
on the values of v and Aq. Figure [3] plots distance be- 
tween the gap of the Leggett mode and the minimum 
energy of the two-particle continuum as a function of A 
for both extended s-wave (T^ = cos k x cos k y ) and stan- 
dard s-wave (Tfc = 1) gaps. 



The principle difference between the two nodeless gaps 
is the range of dopings over which the Leggett mode is ex- 
pected to be observable. In the pure s- wave case, the two- 
particle continuum is given by 2Aq, independent of the 
shape of the Fermi surface. Hence, as seen in Figure [3(b) | 
the dominant effect here is that the gap of the Leggett 
mode scales sublinearly in Ao, and hence becomes ob- 
servable only at large values of the gap. Its separation 
from the two-particle continuum is extremely small at 
the small values of Ao expected to occur near half-filling. 
(Away from half-filling, the gap of the Leggett mode lies 
above the two-particle continuum, as shown in the figure, 
due to the Fermi-surface effects discussed above). In the 
extended s- wave case, however, the minimum of the gap 
also depends on how close the Fermi surface comes to 
the nodes of the gap function. Hence the depe ndenc e on 
filling fraction here is more pronounced (Figure [3(a) | ); for 
all values of Ao the Leggett mode sits definitively below 
the two-particle continuum in the interval 0.45 < v < 0.5. 
Further from half-filling, where the nodes of the extended 
s- wave gap sit closer to the Fermi surface, the mode is 
never visible. Qo is smaller overall in the extended s- 
wave case, compensating for the fact that the minimum 
of the two-particle continuum is de facto smaller than in 
the standard s-wave case. Most notably, for the small 
values of Ao expected near 1/2 -filling, we expect the 
Leggett mode to be below the two-particle continuum in 
the extended s-wave case 
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(a) (b) 

FIG. 3: Distance between the Leggett mode and two-particle continuum for (a) extended s-wave and (b) ordinary s-wave gaps, 
as a function of filling fraction v and gap magnitude Ao in eV. Dark maroon area indicates regions where the Leggett mode 
lies above the two-particle continuum. Blue regions indicate maximum distance from the two-particle continuum. The vertical 
colorbar indicates the energy scale of these differences, as a fraction of the two-particle continuum 2minAo. 



C. Dispersion of the Leggett mode for extended 
s-wave gap 

We now return to the general form of the effective ac- 
tion for the phase degrees of freedom, and analyze the 



structure of its modes at small q. In the absence of the 
Coulomb interaction, the form of the dispersion is effec- 
tively characterized by 



w 



q 2 ^ + c % v ) + n 2 ±J { q \c% v - c| ) + my 



4r 4 



<? 4 cos 2 20 



(13) 



where is the angle in the (k Xl k y ) plane. If vanishes, we retrieve the gapless Goldstone mode, and the gapped 
Leggett mode. From Eq. (12), we find that adding the Coulomb term modifies this according to: 



2 
± 



\J {l 2 (c 2 vv - c 2 H {l - 4N^U(q))) + n 2 } 2 + 



2 Ac% a^l-W^Uiq)) 



cos 2 (20) • 



(14) 



In this case, taking the negative sign for U(q) ~ 
(the unscreened Coulomb interaction in 2D) results in 
uo 2 < 0, indicating that the Goldstone mode has been re- 
placed by a plasma mode. The structure of the Leggett 
mode is, however, largely unchanged by the presence of 
the Coulomb interaction. In particular, we still have 
lim g ^_ uj{q) — VLq. 



Figure 4(a) plots the dispersion relation Q(k) for the 
low-energy modes for several values of Ao, /i. As we have 
kept only terms to quadratic order in we expect 

this to be valid for small q, and have restricted the range 
of the plots accordingly. The important feature to note 
is that the velocity anisotropy, due to the off-diagonal 



terms c^^q 2 — q 2 ), is relatively small and the disper- 
sion is approximately rotationally invariant. Because of 
this, the dispersion relation of the Legge tt mod e is well- 
characterized by the gap VLq (c.f. Sect. Ill A), and the 
velocity v = lim^oo uj(q). This l atter is plotted for the 
extended s-wave gap in Fig. |4(b)[ 



IV. CONCLUSIONS 

We have obtained the fluctuation action for the super- 
conducting phase collective modes of a two-orbital model 
for iron-based superconductors, with particular empha- 
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FIG. 4: Velocity of the dispersion relations Q vs q as a func- 
tion of filling v and gap magnitude Ao for the extended s-wave 
gap A = Ao cos k x cos k y . For small q, the qualitative shape 
of the dispersion (a) (shown at A = 0.1, v — 0.4 ) does not 
depend strongly on either the form of U or the precise choice 
of Ao and v. However, the velocity, plotted as a function 
of Ao and v in (b), is sharply sensitive to the Leggett gap 
and chemical potential. In particular, the magnitude of the 
veolcity increases sharply as the filling fraction decreases. 



sis on the antisymmetric Leggett mode. By fixing the 
parameters of the band structure, our calculation has 
identified the range of doping and superconducting gap 
magnitude over which the undamped Leggett mode ex- 
ists below the two-particle continuum. As the Leggett 
mode's visibility increases with the magnitude of the su- 
perconducting gap, this result also suggests that if the 
bandwidth is narrower, there is a higher possibility of 
observing the undamped Leggett mode. Therefore, a 
strong renormalization of bands could enhance the exis- 
tence of an undamped Leggett mode. Unfortunately, the 
mode and its dispersion are insensitive to the sign of the 
order parameter on Fermi surfaces, and the mode does 
not qualitatively distinguish between the sign-changed s- 
wave and a normal s- wave superconductors. However, we 
find that quantitative characteristics of the mode can in 
principle distinguish between such different pairing sym- 



metries. First, we find that the Leggett mode does lie be- 
low the two-particle continuum, near half filling and suf- 
ficiently deep in the superconducting region, for the ex- 
tended as well as normal s- wave gap. This is distinct from 
the case of nodal gaps, where low-energy quasi-particles 
are always expected to broaden the Leggett mode reso- 
nance. Second, we find that the difference in signatures 
between the two kinds of s-wave pairing symmetry inves- 
tigated here is subtle, but that the extended s- wave gap 
is visible over a narrower range in doping, but further 
below the two-particle continuum over much of its range 
of detectability. This difference comes from the differ- 
ent doping dependence of the two s wave gap functions: 
unlike the normal, sign-unchanging s-wave gap, the s ± 
order parameter will most likely change upon doping as 
the Fermi surfaces become closer to the line of zeroes 
that a sign-changing gap should have in the Brillouin 
zone. This gap variation upon doping is present in both 
strong and weak coupling models^. In this situation, the 
Leggett mode will move from a relatively sharp mode be- 
low the two-particle continuum into a strongly damped 
mode above the two-particle continuum as doping is in- 
creased. This quantitative change can in principle be 
observed in experiments. 



It has been claimed that the Leggett mode has been ob- 
served in MgB 2 by Raman scattering 48 and point-contact 
transport measurement s™, although the energies of the 
Leggett mode measured in the two experiments are differ- 
ent. Ref. 50 found that in the weak-coupling treatment 
of superconductors with an s± gap, however, the A\ g 
Leggett mode does not couple to Raman scattering. The 
analysis carried out here relies heavily on the fact that 
iron-based superconductors are more strongly coupled 
than MgB 2 , and that the superconducting phase is thus 
well-described by considering the orbital, rather than the 
band, basis. This leads to a result which differs from that 
of the weakly-coupled approach in two ways. First, the 
strong-coupling approach suggests that the Leggett mode 
should be observable in Raman spectra. Second, in the 
strong-coupling treatment, the different orbital symme- 
tries should be kept explicitly when determining the rele- 
vant Raman channels. For our model, the Leggett mode 
is caused by an oscillation between the condensates in- 
volving the scattering of a pair of (i^-orbital electrons 
into a pair of d^-orbital electrons. Such a process causes 
a relative density fluctuation, Sn = n xz — n yzi between 
two orbitals, which belongs to the B\ g irreducible repre- 
sentation of the point group (D^h) of the crystals. There- 
fore, the Leggett mode is a B\ g mode in this orbital-based 
model, and should exist in the B\ g channel in Raman 
scattering experiments. (Without the orbital characters, 
the Leggett mode should be a pure A\ g mode, as is the 
case in MgB2^). Thus, observing the Leggett mode in 
channels other than A\ g should provide important evi- 
dence about the orbital structure of condensed pairs in 
the iron-based superconductors. 
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Appendix A: Calculating the phase-only effective action 

We begin by deriving the effective phase-only action for a generic Hamiltonian of the form: 

H = ^2 ^2 C a,r,a( e ( r > r %/3 ~ M<W) C /3,r',a + V a ^(r] r')&£ >r fys,r' + ieU (V - r')p r p r 
a>,(3,r,r' cr=t>i 



(Al) 



Here a, j3 are orbital indices, V is the superconducting interaction (in our model, a spin-spin antiferromagnetic 
interaction decoupled in the Cooper channel to give us cos(k x ) cos(k y ) pairing) , and U is the Coulomb interaction. 
The density is given by p r = ^2 acr ra c a ^ a , and the superconducting bilinear is &£, r 



We may decouple the two 4-fermion interactions by means of two Hubbard- Str at onovich transformations. This 
decoupling gives the action: 



-ie6 r y [ Xr pr> + XrPr'} ~ ^L^a j( r ! r')*/5r' " XrU (r - T>) l \V } 



(A2) 



where <1> is the Hubbard- Str atonovich field associated with the superconducting interaction, and x is associated with 
the Coulomb interaction. 

In computing this effective action, we follow the method of ref. [47lto isolate the action for the phase degrees of 
freedom. That is, taking & ar = A ar e l9a ^ r \ we perform the gauge transformation 



*(r)/2 r 



(A3) 



to absorb all terms involving the phase of the superconducting order parameter into the first term of Eq. ( A2 ) . We 
then have: 



S = E l E C lr,<re- i9a(r)/2 [(d/dr " M)<W + ^(r, ^)] <*,r> e"^'" 2 - S r y 
r,r f v <y 

?5 r , r > [ Xr pr> + XrPr'] ~ $L^cJ( r ; r ')*/3r' " XrU(r - r'^Xr'j 



a.,f3,r,r' 

—ied 



(A4) 



It is convenient to re-express the kinetic terms as: 

Ea, W { 4,r,t e "^ (r)/2 " M)<W + e«Ar, r')] c^ A e i6 ^' 2 



-C a ,r,ie- ie ^' )/2 [(-d/d T ~ p)Sa,0 + e^V)] C^e"*™ 2 

+5 a Ar,r>e- i(>air)/2 \{d/d T - /i) + W*(r', r)] e i0 ^' 2 ) 



(A5) 



The first two terms can now be combi ned with the rest of the fermionic action can be expressed in matrix form in 
the BCS basis. The final line of Eq. (A5) gives a separate contribution to the action, which can be expressed, to 
quadratic order in 6,d r 6: 



d ri d r , 



)} 



(A6) 



The first term is a total derivative and will not contrubute to the dynamics of the phase-only effective action. The 
second term in any case vanishes, as v is odd over the Brillouin zone. Hence only the last term appears in the effective 
action. 
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As we are interested in the dynamics of the phase degrees of freedom, we replace A a?r by its mean-field value. For 
the time being, we drop the Coulomb terms by setting U = 0, x = 0; we will discuss these further in Sect. [B] We 
further consider only the slowly varying phase fluctuations. This allows us to expand the exponentials of the fermionic 
terms in (A4) using: 



e -i(Mr)-Mr')) ~ i _ * (0 a(r ) _ 6a {r')) - l(6» a (r) - e a {r') f 



e 2 



(9l(r)-6 2 (r>)) 



e -i<po | 



1 - |(0i(r) - 9 2 (r') - <p ) - \{6 a {r) - 6 a {r') - J 



(A7) 



where we have explicitly separated out the possible background expectation value of the phase difference ipo between 
the gaps. In practice cpo = is set by the mean-field equations. 

Defining = (c[ r ^c\ iT ^c! 1 r ^,C2,r,i), we may now express the first two terms in Eq. (A5), after Fourier trans- 
forming, as 



s f 



du\du 2 



d 2 ki d 2 k 2 ,+ 

oh < 



47T 2 47T 2 



with 



G 1 (k u k 2 ,iuj u iuj 2 ) = G 1 (k 1 ,oj 1 )5 qj0 Sn j o + k 2 , iu>i, iv 2 ) 
where here q = k\ — k 2l Q = u)\ — oj 2 . We have: 



I iuj + e x Ai 
Ai iuj — e x 

xy 

\ o 



^xy 









^xy 

A 2 iu — e 



(A8) 



(A9) 



(A10) 



y 



where e a = e a (k) is the kinetic energy in the orbital basis, and A a = A a (k) is the momentum-dependend supercon- 
ducting gap in each orbital. The second part of (A9) is given by: 



q, iuo, id) 



n f e 1(Jz o 

2 2 a z 



1 



i_ ( 1 S q (e x )a o fatxyifa) - #ie^(/c 2 )] <r 

2 V [®i e xy(ki) ~ 2 e xy (k 2 )} a 2 S q (e y )a 

r(2) 



^ (#i(/c 3 ,^ 3 ) 2 (g-& 3 ,i(Q-cj 3 )) 



/c 3 ,2u;3 











i(q - k 3 ,i(Q - uj 3 )) 
2 (q-k 3 ,i(tt-uj 3 )) 



B{k\,ibJ\, k 2 , iuo 2 , k 3 .iu 3 )a z 

C(fci, zgji, fc 2 , ^2, k 3 .iuj 3 )cr z 



(All) 



where 0q, = 6 a (q,iQ), and we have defined the discrete derivatives: 



<Va/3 = e a/ 3(fci) - e a/3 (/c 2 )) 
4 2)e a/3 = e a)9 (fei) - e a(3 (k 2 + k 3 ) - e ai9 (fci - fc 3 ) + e Q , /3 (/c 2 ) 



(A12) 



The off-diagonal terms quadratic in the phases are: 



B(ki,iu)i, k 2 ,iuj 2l k 3 ,iuj 3 ) = (0 2 {k 3 ,iuj 3 )0 2 (q - k 3 ,i(Q - uJ 3 ))e xy (k 1 ) - 2 (k 3 ,iuj 3 )0 1 (q - k 3l i(Vt - uj 3 ))e xy (k 2 + k 3 ) 

-0 1 (ks,ius)6 2 (q - k 3 ,i(Sl - uJ 3 ))e xy (ki - k 3 ) + 0i(fc 3 ,io; 3 )0i(g - k 3 ,i(tt - u 3 ))e xy (k 2 )) 

C(ki,iu)i, k 2 ,iu 2 , k 3 ,iuj 3 ) = (0 1 (k 3 ,iu 3 )0 1 (q - k 3 ,i(Q - u 3 ))e 21 (k 1 ) - 0i(fe 3 , zcj 3 )0 2 (g - fc 3 ,i(fi - u 3 ))e xy (k 2 + fc 3 ) 

-0 2 (fe 3 , wj 3 )0i(g - ks,i(Q - uj 3 ))e xy (k 1 - k 3 ) + 2 (k 3 ,iu 3 )6 2 (q - k 3: i(tt - u 3 ))e xy {k 2 )) 

(A13) 

Thus in our treatment, the block diagonal terms involve only discrete differences of the band energies, which will 
become derivatives when the momentum of the phase variables is small. The off-diagonal terms contribute, as well as 
such differences, a term which is finite at q = (or k\ = k 2 ). Hence the gap of the Leggett mode is, in the absence of 
inter-orbital pairing, generated by the kinetic mixing between the two orbitals. 



12 



To obtain the effective action, we integrate out the fermions in Eq. (A8). In practice, we must evaiuate the resuit 
perturbativeiy in E. Specifically, we have: 



?e// = 5 M F-Trln(l-G S) 

w Smf + Tr (G S) + \tt (G £G S) 



(A14) 



where Smf is the mean-field action, from which we self-consistently determ ine the values of Ai, A2. Here we will 
evaluate the low-energy, long- wavelength limit of the effective action (A14) by keeping terms to quadratic order in 
q, and a (q, Q). 



1. Evaluating TrGoT, and TrG^G^ 

For reference, here we give a more detailed account of the calculation in Sect. [TTJ 
The separate expressions for the two traces are: 



Tr(G £) 



-/ 



d 2 k 



(ip(q,in)ip(-q,-iQ)) 



(2tt)2 I . 



1 



E+(E+-e y ) 

? (A) 



E-(E--ey) 



(E+-E- ) E K _> (E+-E- ) 



,(12) 

L ij (E+-E-) 



? (A) 



{E + -EJ)Ef ) E. 



€ *y ( ei(a) 77- 



(2) 



(12)_e £L 
x _ E + (E + -e x ) 

e[ a \e + -e_) 



p(A) 



p(A) 



r (A) 



/ 



(A15) 



?2(-(7) 



where m. 



_ Q k € Qfc. •> anc ^ we define (/9(g) = 6^1 (<^) — 62(0)- We have dropped the linear term in ^, because it is a total 
derivative and hence should not contribute to the action. Here all e, E, and A are evaluated at the momentum k to 
be integrated over. Note that we have also included the quadratic terms in the last line of Eq. (A5). 
Evaluating Tr(GoXGoS) gives: 



1 



Tr(GEGE) 



n 2 



/ 



( <f%) <p(q) 



N 6 







N, 



LpLp 



<p(-q) 



d 2 k 

WT 2 



• e 2 -E (A) E (A) 



■ E+E- 



(a) e (a ){e (a ) + e (a: 

A 2 + E + E_ \ 



1 



8 (4 A) + E W )(E + - E_) 2 y (E^E™ 
( ~ ( e ly + ~ e y?) v\ xy) vf y) - e 2 xy v i4> v j4> 



( HqMq) ) 



(Vi(/>Vj<p + Vj^Vi^) 



f- (v i4> v j(p + v j<f> v i(p ) 



+ 



^xy (tec 





A(f)ip A^ptp 



e v ) ( (xy) 
_y_L \ v \ yJ v 

\\( H-q) 



W , (xy) (4>)\ 
3 ~ 3 i I 



4 



(e 2 xv -W*-ty) 2 )4 Xy) 4 Xy) 



e_A (v\ XV) V i 



3 



J 1 



(A16) 



where = d(e x + e y )/dki, v^i = 9(e x — e y )/dki. Here A a £ are terms which come from expanding traces involving B 

in Eq. |A13| to quadratic order in q. 

Combining the mass terms from Eqs. ( |A15| ) and (A16) gives the total mass term: 



M 



d 2 k 
(2^) 2 



4A 2 e 2 

^xy 



e{ a) e (_a) {e (a) 



^^xy 



d 2 k 



4A^e 



2^2 



xy 



(27r)2 £ ;{A) f ,(A) (£7 {A) +£ ,{A) ) 



(A17) 
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where the second equality holds because in practice e xy averages to over the Brillouin zone. 
In simplified form, the momentum-dependent terms are: 



N c 2 

1 v <£>wj 



-II 



d 2 k 

(2tt)2 



(E+ - E. 



E. 



E. 



(A) 



(xx) (xy) . (yy) (xx) (yy) 

m\j U x -2m\ j y >e xy +m) yy >e y m\j >e x - m\ yy >e y 

(xx) (yy) (xx) (xy) (yy) 

'r mi" <„ n>, :! 'e x + 2m\ j y! e xy + rnf j y 'e y 



{Ef ] + E^'){E+ - E. 



p(A 



(A) 771(A) 



(Ef'E. 







(xy) v W 



e xy {e x -e v ) 
4 



{el y -\{e x -e y f)v^v^- 



(xy) ((f)) . (xy) (0)' 
J J 1 



A 

A^ A 



(pep 



Appendix B: Effective action with Coulomb terms 



Inc luding terms generated by the Coulomb repulsion modifies the interaction term E of the full fermion propagator 
(All) according tcH3 : 



E = r 3 <g> l(z— zex) + r 3 r 3 i— - + E A 



kin 



(Bl) 



where E^ n involves only spatial derivatives of the phases 0\ and #2- Here x is the Hubbard-Stratonovich field 
associated with the Coulomb interaction. The form of the coupling for \ to fermions can be deduced from gauge 
invariance: the phase 6i is obviously a gauge-dependent quantity, and the gauge- invariant degrees of freedom are 
the combinatio ns d T 0/2 — e\ — eAo and X/0/2 — e/cA^A Hence the effective action for x is the same as that for 
d T 0/(2e). Eq. (Bl) indicates that x couples in all cases like the time derivative of the symmetric component of the 
phase fluctuations. 

To obtain the full effective action for the phase only modes in the presence of Coulomb interactions, we first 
integrate out the fermions, giving an effective action for the 3 Hubbard-Stratonovich fields </>, (/?, x- There are two 
relevant contributions: from TrGE, we obtain: 



ie X (q)Tr [(r 3 l)G k . 



-iex(q)(pk- q ) 



(B2) 



which cancels the first-order term in x i n the effective action (A2). 

From TrGEGE, we obtain contributions whose coefficients are the same as the contributions from the time deriva- 
tives of (j). In particular, as the coefficients of the cross-terms in q, Q from traces GEGE all vanish, the couplings 
between x an d </>, ip depend only on Q. Hence the effective action for the fields (/>, tp, x nas the form: 



dQd 2 q ( </) 



n 2 - 



N vv [9? 



Ml - qD 




(B3) 



where the coefficients A/", c are given in Eqs (|5|, ( |A15 ), and ( A16). The dispersion relation is given by finding the values 
of g, Q at which A4 is singular. Depending on the values of the parameters, M. may have one or two modes which 
are finite as q —> 0. One of these is the gapped Leggett mode; the other is a sound-like mode (the Carlson-Goldman 
mode) which we find to be absent at T = 0, consistent with Ref. [47j The third mode is, of course, the plasma mode, 
which does not appear in the low-energy spectrum. To study only the phase modes, we may equivalently integrate 
out x an d 4> to obtain Eq. (11 ). 
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